Localization in an Inhomogeneous Quantum Wire 
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We study interaction-induced localization of electrons in an inhomogeneous quasi-one-dimensional 
system — a wire with two regions, one at low density and the other high. Quantum Monte Carlo 
techniques are used to treat the strong Coulomb interactions in the low density region, where 
localization of electrons occurs. The nature of the transition from high to low density depends on 
the density gradient — if it is steep, a barrier develops between the two regions, causing Coulomb 
blockade effects. Ferromagnetic spin polarization does not appear for any parameters studied. The 
picture emerging here is in good agreement with measurements of tunneling between two wires. 



With the rapid development of nanotechnology over 
the last decade, experiments have been able to probe 
strong interaction phenomena in reduced dimensional- 
ity systems such as quantum dots, wires, and point con- 
tacts. Of particular interest are systems in which the 
electron density is inhomogeneous. In the low density 
region of such a system the interaction energy is com- 
parable to the kinetic energy and novel effects occur 
such as the "0.7 structure" in quantum point contacts 
or Coulomb blockade effects accompanying localization 
in a one-dimensional (ID) wire. We perform quantum 
Monte Carlo calculations of an inhomogeneous, quasi-lD 
electron system in order to address such effects. 

The so-called "0.7 structure" [l|] remains poorly under- 
stood: in a quasi-lD electron gas - a wire, constriction, 
or quantum point contact (QPC) - decreasing the density 
causes the conductance G to decrease in integer multiples 
of 2e^ /h (one for each transverse mode), except for an ex- 
tra plateau or shoulder at G « 0.7(2e^/h) as the lowest 
mode is depopulated [2, Isl, \4, la, la, lZ|, l8| . Proposed theo- 
retical explanations have been mainly based on three ap- 
proaches: formation of a bound state leading to a Kondo 
effect [9|, llOj , spontaneous s pin polarization of the low 
density electroiis M , M M, III, llll, ll2] , and formation of a 
Wigner crystal [13|, llJ] . An open question is whether the 
critical features underlying each of these approaches is 
present in an inhomogeneous quasi-lD system — whether 
a localized state with Kondo-like correlations, spin po- 
larization, or a Wigner crystal occurs. 

The formation of a Wigner crystal was investigated di- 
rectly using_tunneling spectroscopy into a quantum wire 
ISl . lid ll7l . Il8( . Clear evidence of localized electrons was 
found, accompanied by unexpected single electron phe- 
nomena. 

The general problem of a transition from a liquid to 
a localized crystal-like phase remains a subject of funda- 
mental research in a va riety o f bulk and nanoscale sys- 
tems [UliSlIl HI, H, 



Thus, an approach from 
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of interaction-induced liquid to crystal transitions. 

Previous electronic structure calculations investigating 
inhomogeneous ID systems have been based on mean 
field approximations. While some density functional cal- 
culations in the local spin density approximation (LSDA) 
were used to support the spontaneous spin polarization 
other LSDA and Hartree-Fock (HF) 
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scenario 



calculations 



27l . |28| . |29| confirmed the existence of 



a quasi-lD point of view is valuable not only for under- 
standing the 0.7 anomaly and the tunneling experiments, 
but also for bringing a new way of looking at the physics 



quasi-bound states, which may lead to Kondo-like physics 
in a short QPC 1J, ISj. A very recent HF calculation in 
a long constriction in a weak magnetic field predicted an 
antiferromagnetic to ferromagnetic transition [30|. De- 
spite considerable interest, no many-body calculation be- 
yond mean-field has, as far as we know, been performed 
for an inhomogeneous ID system. 

Here we use variational and diffusion quantum Monte 
Carlo (QMC) techniques to investigate correlation and 
localization in a zero temperature, inhomogeneous, quasi- 
lD electron system. We observe interaction-induced lo- 
calization in the low density region. Our focus is on 
the transition between the high and low density regions, 
which can be either smooth or abrupt and may involve 
the formation of a barrier, depending on the smoothness 
of the external potential. 

To study interacting electrons in an inhomogeneous 
quasi- ID system, we consider a narrow two-dimensional 
quantum ring with a constriction (point contact): 

N N N 

i i i<j ■^ 

+ Vg{tanh[sie, + eo)]-tanh[s{e,-eo)]} (1) 

where effective atomic units are used — the length scale 
is Oq = h^e/rit^e^ and the energy scale is the effective 
Hartree H* = e^/eag. For comparison with the tunnel- 
ing experiments, the values for GaAs are a^ — 9.8 nm and 
H* = 11.9 meV. The parabolicity oj controls the width of 
the ring, vq is the radius, and Vg is the gate voltage that 
controls the electronic density in the low density constric- 
tion. The sharpness and length of the gate potential are 
tuned through the parameters s and ^o- In this work, uj 
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FIG. 1; (Color online) Two-dimensional ground state density 
p for gate potentials of different shape [all with Vg = 0.8 (in 
units of H*) and TV = 31]. Potentials A, B, and C have 
s = 15, 4, and 2, respectively. D is obtained using a gaussian 
with an angular width of 1.2 radians. Potential B is used in 
most of the paper. 



and ro are set to 0.6 and 25; 6*0 is 1.5 radians. Our po- 
tential yields a low density region of length 75 (0.73 ^m 
for GaAs). iV = 30 or 31, generating a single-mode ID 
electron density of n ~20 ^m~^ when V^ = 0. The elec- 
tron gas parameter, Ts = l/2naQ, is 2.6 in this case. For 



previous work on (uniform) quantum rings, see Ref.l31. 

Variational Monte Carlo (VMC) is the first step of our 
numerical approach. We use recently developed energy 
optimization methods |32l . |33[ to optimize a Jastrow- 
Slater type trial wave function *t(R.) = J(R)L'(R) [si ]. 
To build the Slater determinant D(R), we considered 
three qualitatively different types of orbitals: x(^) e"*"*"^ 
with x{i") fixed, LSDA orbitals, and floating gaussians 
[22j . Most of the results presented here were obtained 
using floating gaussians as they provide a better descrip- 
tion in the strongly localized regime - a comparison is 
given below. After optimizing the variational parameters 
(Jastrow parameters as well as, separately, the positions 
and radial/ angular widths of the floating gaussians), we 
then perform a diffusion Monte Carlo (DMC) calculation 
to project the trial wave function onto the fixed-node ap- 
proximation of the true ground state [35| . The fixed-node 
DMC energy is an upper bound on the true energy and 
depends only on the nodes of '^t obtained from VMC. 

The density plots of Figs.[l] and [5] give an overview of 
different scenarios that can occur in an inhomogeneous 




FIG. 2: (Color online) For a short constriction, the two- 
dimensional ground state density, showing a single localized 
electron. [Vg = 0.5 (in units oi H"), Oq = 0.7, and N = 30.] 



quantum wire depending on the gate potential landscape. 
Potential A, which is close to a square barrier, gives rise 
to three interesting phenomena: (i) The ripples in the 
high density part of the ring are Friedel oscillations (1 
maximum per 2 electrons). This is a signature of weak 
interactions - the electrons are in a liquid-like state, (ii) 
The modulation in the low density part of the wire shows 
4 electrons that are individually localized, (iii) There is a 
large gap separating the liquid and crystal phases, caus- 
ing the low density region to be in the Coulomb blockade 
regime. This effect has been observed experimentally. 



!rved experimental 



but the origin was not understood 

For a smoother potential step, the size of the gap de- 
creases (potential B in Fig.[T]) and eventually disappears 
(potential C). Finally, for the gaussian-shaped potential 
D, localization is very weak. 

As the length of the low density region is reduced, the 
number of localized electrons decreases. In the extreme 
case, the low density region becomes like the saddle po- 
tential of a quantum point contact. Fig. [2] shows that for 
a short constriction, a single electron is localized in the 
low density region, with substantial barriers to the high 
density "leads", as in LSDA calculations [1,0, [S^]. The 
spin of the localized electron fluctuates, as in the "Kondo 
scenario" 0,11, [i3| for the 0.7 structure. 

The results in Figs.[l]and[2]show that (1) it is possible 
to have localized electrons in a low density region dis- 
tinctly separated from the liquid leads and (2) an abrupt 
potential barrier with a flat plateau enhances both local- 
ization and the gap between liquid and crystal regions. 

Fig. [3] presents how localization develops as the gate 
potential increases. For the rest of the paper we fo- 
cus on potential B {s — A), for which the width of the 
potential riser is ~ 60 nm. When Vg = 0.4, the den- 
sity modulation indicates that localization of electrons 
is already beginning. The average density at the con- 
striction is ~ 15 iim.~^ {rg ~ 3.4), which is close to 
but lower than the experimentally estimated critical den- 



sity n* 



20 /im for localization [15[. For small Vg 



{Vg < 0.2), no oscillations occur. On the other hand, 
as Vg increases, the electronic density decreases. The 
liquid/crystal gap becomes clearly visible at Vg = 0.7. 
At that point 5 individually localized electrons are well 




FIG. 3: Electron density as the low density region is depleted. 
Vg (in units of H*) is varied for potential shape B (s = 4, N = 
31); each curve is shifted vertically by 20 fj,ra~^ with respect 
to the one below. For Vg — 0.4, 8 electrons are localized. As 
Vg increases, localization becomes stronger and a gap forms 



at the crystal-liquid boundary- 
blockade regime. 
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separated from the bulk {ts ~ 5.7). Further increases 
in Vg cause the number of electrons in the constriction 
to decrease abruptly at certain values while the liquid- 
crystal gap becomes even stronger. The physical picture 
that emerges from our QMC calculations is remarkably 
similar to that in the momentum resolved tunneling ex- 
periments [il,[il,[i3,[ii. 

To investigate the formation of a barrier between the 
liquid and crystal regions, consider a model calculation: 
suppose the gate voltage has depleted the low density 
region to the point that no electrons reside there. Clearly, 
an additional electron in the low density region would 
feel a barrier caused by the electrostatic repulsion of the 
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FIG. 4: (Color online) The total effective potential (electro- 
static plus external) of bulk electrons in a square well [poten- 
tial is 1.8 for a; > and x < —150, otherwise (other param- 
eters given in [36])]. Solid line is from a quantum mechanical 
single-electron calculation (Inset: density, with Friedel oscil- 
lations comparable to those seen in Fig. [T] for potential A). 
Dashed line is obtained assuming a constant density distri- 
bution. Clearly, an electron at x > would feel a potential 
barrier separating it from the bulk electrons. 



FIG. 5: (a) Difference of total energy when electrons in the 
constriction are fully polarized or antiferromagnetically or- 
dered, as a function of Vg {N = 31). The antiferromagnet is 
lower in energy for Vg < 0.6. (b) For a homogeneous ring with 
30 electrons, energy difference per electron between fully po- 
larized and antiferromagnetic states. The antiferromagnetic 
state is always favored for the parameters we have studied. 
{Vg = 0; density varied by changing ro.) 



electrons in the leads. For example, Fig. 0] shows for two 
situations the effective potential due to the repulsion of 
the bulk electrons plus a square barrier gate potential 
(finite step at x = 0). Our QMC results suggest that a 
barrier is similarly formed even when the density in the 
constriction is not so depleted, though the discreteness of 
charge, quantum interference, and correlation probably 
play a role in enhancing the effect. 

The spin structure of the electrons in the low density 
region is an irnportant physical pro perty which has been 



controversial 
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In one 



dimension, the ground state cannot be ferromagnetically 
polarized (Lieb-Mattis theorem [37|), but the transverse 
degree of freedom may change the situation. With our 
method, a rigorous treatment of the spin is technically 
difficult: building eigenfunctions of S^ out of floating 
gaussian orbitals for a large number of electrons with 
S <C N/2 requires a huge number of determinants. We 
can solve this problem by doubly occupying extended or- 
bitals (i.e., x(^) e^™* or LDA orbitals if a self-consistent 
solution can be found), but this is significantly less ac- 
curate than floating gaussians. Thus, here we present 
results using Sz conserving floating gaussian trial wave 
functions. 

First, we find the energy difference for a homogeneous 
ring (i.e., Vg = 0) between fully polarized and antiferro- 
magnetic states. Fig.ISJb) shows that the ferromagnetic 
state is always higher in energy and that the difference 
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FIG. 6; (Color online) Comparison of ground state densities 
obtained using LSDA, QMC with LSDA orbitals, and QMC 
with gaussian orbitals at (a) Vg — 0.7 and (b) Vg = 0.87. 



decays exponentially in -^/rj, as expected |13l.ll4|. 

In order to compare to experiment [l5|, we also study 
the spin state in the inhomogeneous ring. Fig.[5{a) shows 
the energy difference between ferromagnetically and an- 
tiferromagnetically arranged electron spins in the con- 
striction. (In the high density "leads" , the spins have no 
definite arrangement because of strong fluctuations com- 
ing from weak correlations.) At Vg — 0.5, the low density 
part is clearly not fully spin polarized as the spin gap is 
as large as ~ 3 K. However, as Vg increases, both the 
number of electrons and their overlap decrease rapidly, 
resulting in a much smaller energy gap. For Vg > 0.6, 
the energy gap is not resolved due to statistical error of 
about 0.5 K. This is consistent with Steinberg et aL's ex- 
periments performed at temperatures down to ~ 0.3 K 
[15|: it was found that for iV < 6 the observed state is a 
mixture of ground and thermally excited spin states. 

We close by comparing the results obtained with dif- 
ferent methods and trial wave functions. Fig. [6] shows 
the densities obtained from LSDA and from two QMC 
calculations, one using LSDA orbitals in the trial wave- 
function and the other using floating gaussian orbitals. 
In Fig.inja) there is excellent agreement between the two 
QMC calculations whereas the LSDA density is not as 
localized. Fig.[S]Jb) shows that the inaccuracies of LSDA 
become more visible at higher Vg (lower density), as ex- 
pected. The gap that forms between the localized and 
liquid regions is bigger in QMC than in LSDA, indicating 
a correlation contribution to the gap. We also performed 
QMC calculations using x{i") e^*"^ orbitals (hence fixing 
the total spin), but the density is substantially different 
and the energy considerably worse. Finally, although the 
fixed-node DMC energies obtained from LSDA and gaus- 
sian orbitals are very close [51.7915(6) and 51.7862(2) re- 
spectively at Vg = 0.87], the floating gaussian based trial 
wave functions yield significantly reduced fluctuations of 
the local energy (standard deviation of 0.45 compared to 



0.59). Thus floating gaussian trial wave functions pro- 
vide a better physical description of the system in the 
Wigner crystallized regime. 
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